Impact of nuclear "pasta" on neutrino transport in collapsing stellar cores 
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Nuclear "pasta", nonspherical nuclei in dense matter, is predicted to occur in collapsing supernova 
cores. We show how pasta phases affect the neutrino transport cross section via weak neutral current 
using several nuclear models. This is the first calculation of the neutrino opacity of the phases with 
rod-like and slab-like nuclei taking account of finite temperature effects, which are well described 
by the quantum molecular dynamics. We also show that pasta phases can occupy 10-20 % of the 
mass of supernova cores in the later stage of the collapse. 

PACS numbers: 26.50.+X, 25.30.Pt, 97.60.Bw, 13.15.+g 



Interactions between neutrinos and matter are one of 
the key elements for understanding the mechanism of 
core-collapse supernovae. In the collapse phase, neutrino- 
matter interactions play a crucial role for explosion: the 
efficiency of the trapping of neutrinos inside the collaps- 
ing core governs its lepton fraction Yt. (the ratio of the 
number density of leptons to that of baryons) and hence 
the mass of the inner core at bounce. The latter quan- 
tity determines the initial position and the strength of 
the shock wave [l| ■ 

In the neutrino-trapping stage, interactions between 
neutrinos and matter are dominated by coherent scatter- 
ing of neutrinos from heavy nuclei via the weak neutral 
current, which we focus on in this article. In the case 
of scattering of neutrinos by an isolated nucleus, neutri- 
nos are coherently scattered by nucleons in the nucleus 
provided the wavelength of neutrinos is much larger than 
the radius of the nucleus r^. Therefore, roughly speak- 
ing, the cross section per nucleon is enhanced by a factor 
of A 0, H|, where A is the mass number of the nucleus, 
due to the fact that the scattering is coherent. How- 
ever, when neutrinos are scattered by matter in the core, 
not by an isolated nucleus, more than one nucleus would 
contribute to the scattering amplitude and the coherence 
can be disturbed if the internuclear distance is compa- 
rable to or smaller than the wavelength of neutrinos [J] . 
This correction is due to the spatial correlation among 
nuclei. Especially, at higher densities, nuclei form a lat- 
tice because of the strong Coulomb forces among nuclei 
and thus the correlation among them is very different 
from that of the liquid phase at lower densities. 

More interestingly, the nuclear shape itself would 
change from sphere to rod or slab, etc., at subnuclear den- 



sities in the collapsing supernova core [f| @, 0, H, H, EH • 
These nonspherical nuclei are referred to as "nuclear 
pasta" . Transitions to pasta phases also change the cor- 
relation among nuclei and hence the neutrino opacity. 

At the present time, the effects of nuclear pasta on the 
neutrino opacity have yet to be understood completely. 
One of the important questions to be clarified is how the 
neutrino opacity for the conventional models of dense 
matter used in supernova simulations is modified by tak- 
ing account of the pasta phases. Horowitz et al. have 
studied a related problem using the framework of quan- 
tum molecular dynamics (QMD) pT| , but they have not 
reproduced the phases with rod-like and slab-like nuclei 
and those with rod-like and spherical bubbles probably 
due to the use of a short screening length (10 fm) com- 
pared to the internuclear distance (see a discussion in 
Ref. [13] )■ Therefore the above question is still open. 

We calculate the cross section for elastic neutrino scat- 
tering from pasta phases via weak neutral current [HI]. 
After averaging over the spin states of nucleons, result- 
ing expression of the neutrino cross section per nucleon 
to the first order is 
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where N is the number of nucleons in the system, Gf is 
the Fermi coupling constant, E v is the energy of neutri- 
nos, 6 is the scattering angle, cl™" 1 (c^) is the vector cou- 
pling constant of neutrons (protons) to the weak neutral 

current (cl"' = -1/2, cl p) = 1/2 - 2 sin 2 9yv, sin 2 9 W = 
0.23), hq is the momentum transfer, and S(q) is defined 
as follows: 
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In Eq. x p (x n ) is the proton (neutron) fraction of the 
system, Sab(<1) with A, B — n or p is the partial static 
structure factor of nucleons of species A and B, 



Sab (q) = 1 + p / d\ e iq r { 5AB (r) - 1}, 



(3) 



where p is the average number density of nucleons in the 
system and gAB is the two-body distribution function of 
species A and B defined as follows, 
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In Eq. ([3]), (• • • ) denotes the statistical average, tpA (V's) 
is the field operator of particle A (B), pa [pb) is number 
density of A (B), and V is the volume of the system. If 
we ignore the contribution of the axial vector current, the 
cross section for the scattering of neutrinos by a single 
neutron is given by 
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Comparing Eqs. fl} and ([5|), we see that S(q) is the am- 
plification factor of the cross section per nucleon due to 
the structure of nuclear matter. We normalize Eq. @ 

by c« because neutrons dominantly contribute to the 
vector current. 

In the field of simulation of supernovae, researchers 
often use transport cross section defined as 
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which means the efficiency of the momentum transport 
from neutrinos to matter. Note that the cross section de- 
pends on the relative direction between the momentum 
of neutrinos and the symmetry axis of the pasta phases 
as shown later in Eq. (|10p . However, collapsing cores 
are polycrystalline; the symmetry axes of grains are ran- 
domly oriented and the size of the grains would be much 
smaller than the mean free path of neutrinos. Thus the 
neutrino opacity of supernova cores would be well char- 
acterized by the angle-average of S(q) (an average over 
the scattering angle 9), 
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where Qi and Qf are integrated over initial and final di- 
rections of neutrinos. The above quantity is related to 
the transport cross section a t per nucleon as 



a t = (S{E u ))al 



where a® is the total transport cross section of one neu- 
tron via the vector current, 
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Thus (S(E V )) is an amplification factor of the total trans- 
port cross section. 

When nuclei form a crystalline lattice, the partial static 
structure factor 5 J 4s(q), within the Wigner-Seitz (WS) 
approximation, can be written as 
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where G is the reciprocal lattice vector of each structure, 
xa {xb) is the fraction of particles of species A (B) in 
the system, Vws is the volume of the WS cell, Fa {Fb) 
is the form factor of a single nucleus for species A (B) 
defined as 
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Here p^ is the difference of densities of species A inside 
and outside nuclei, and the integration J v is over the vol- 
ume V n of one nucleus. In the following, we use the above 
expressions to calculate (S(E l/ )) for models given by our- 
selves [H], EH and by Lattimer and Swesty jiH | in which 
lattice structure is treated within the WS approximation. 
In the present calculation, however, we take account of 
the lattice structure explicitly: we assume a BCC lattice 
of spherical nuclei and bubbles, a triangular lattice of 
rod-like nuclei and bubbles, and a one-dimensional lat- 
tice of slab-like nuclei. These assumptions are justified 
by the fact that the Coulomb energy for the assumed 
lattice structure is well reproduced by the WS approxi- 
mation [l7[. In the following numerical calculations, the 
summation about G in Eq. (|10p is done for the six small- 
est reciprocal lattice vectors because the contribution of 
higher orders is negligible in the present situation, where 
the nuclear volume fraction is large (the contributions of 
the fifth and sixth orders are already 10~ 3 smaller than 
that of the first order) . 

To see the effects of existence of nuclear pasta clearly, 
we compare (S(E V )) for two cases: (1) we take account 
of the pasta phases, i.e., the phases with rod- like nuclei, 
slab-like nuclei, and rod-like bubbles, in addition to the 
phases with spherical nuclei and bubbles; (2) we only take 
account of the phases with spherical nuclei and bubbles 
(here referred to as the no pasta case). For case (1), we 
use nuclear data of supernova matter at zero temperature 
and at Yl — 0.3 obtained in Refs. [3] and [Hj]. In these 
references we used a liquid drop model, which is based 
on that developed by Baym, Bethe, and Pethick (BBP) 
[l8l |. Here we use results for standard values of param- 
eters for the surface tension and the chemical potential 
in this model. Case (2) is newly calculated for the same 
parameters using the same model and taking account of 
only spherical nuclei and bubbles. In Fig. Q] we show the 
results of these two cases. We can see the common ten- 
dencies that (S(E V )) is zero below E v ~ 25 MeV, rapidly 
increases at 25-30 MeV, has a peak around 30-40 MeV, 
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FIG. 1: (Color online) (S(E„)) at Y L = 0.3 calculated for a 
liquid drop model based on BBP with a typical parameter set 
given by Ref. 0. The abbreviations SP, C, S, CH, and SH 
stand for phases with spherical nuclei, cylindrical nuclei, slab- 
like nuclei, cylindrical holes and spherical holes, respectively. 
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FIG. 2: (Color online) (S(E„)) at x p = 0.3 calculated for EOS 
by Lattimer and Swesty (LS-EOS) 0. In the left panel, the 
temperature is fixed at 1 MeV. In the right panel, the number 
density of nucleons is fixed at 0.0165 fm -3 . 



falls relatively gently and approaches x n + (c« ' /cv) 2 x p 
above 100 MeV Peak height monotonically de- 

creases with increasing density. This is simply because 
the nucleon density distribution becomes more uniform. 
Striking differences between the results for the cases with 
and without the pasta phases are as follows. The peak 
height decreases and the peak energy increases when one 
takes account of the pasta phases. The former is caused 
by the disappearance of the coherence along the axis of 
cylindrical nuclei and in the plane of slab-like nuclei. The 
latter is due to smaller values of the lattice constant for 
the phases with rod-like and slab-like nuclei compared 
to those for the phases with spherical nuclei and bub- 
bles; this is a general property observed in a number of 
previous calculations 0, EE E3, EE EI . 

In Fig. [2] we show (S(E„)) for the equation of state 
(EOS) by Lattimer and Swesty (LS-EOS) [Hj] at a pro- 
ton fraction x p = 0.3 to see effects of finite tempera- 
tures within the liquid drop approach 22]. The left panel 
shows the density dependence and the right one the tem- 
perature dependence. We note that, in this calculation, 
nuclei are assumed to be spherical and the other struc- 



tures including the spherical bubbles are not taken into 
account. In the present calculation, when we take an av- 
erage of 5,4s (q) over the scattering angle, we multiply 
5y4s(q) by the Debye- Waller factor given by Ref. [231 ] 
and use the relation q 2 = 2E 2 (1 — cos &). The right 
panel shows that (S(E V )) falls monotonically with in- 
creasing temperature. This is due to thermal fluctuations 
of the position of nuclei from equilibrium lattice points 
described by the Debye- Waller factor and the smoothed 
nucleon density profile at higher temperatures. Although 
these two kinds of finite temperature effects can be (even 
partially) incorporated within the liquid drop model such 
as LS-EOS, there is another kind of finite temperature ef- 
fect that is hard to be described by the liquid drop model: 
disorder. Figure [5] exhibits a steep decline of the peak on 
the low energy side. Disorder will broaden the peak of 
(S(E V )) and will produce a spectral tail on the low energy 
side since the disorder of lattice contributes to (S(E V )) at 
wavelengths necessarily longer than the lattice constant. 
The QMD can describe the finite temperature effects in- 
cluding disorder beyond the limit of liquid drop models. 
Next, we show (S(E U )) calculated for QMD simulations 
as more realistic results. 

Here we briefly summarize QMD models and com- 
putation techniques. We use nuclear forces given by 
Maruyama et al. (24| (model 1) and by Chikazumi et 
al. [IE and Kido et al. [IE (model 2) with medium 
EOS parameter sets \ 
ous works H EE HE E 



. Model 1 is used in our previ- 
Model 2 is newly adopted for 
studying the dependence of the phase diagram on nuclear 
forces within the framework of QMD [^E ■ Parameter sets 
for both the models are determined by fitting the satura- 
tion properties of nuclear matter, the properties of finite 
nuclei in the ground state, and the energy dependence 
of the optical potential. In the Hamiltonian of model 2, 
there is an extra term depending on the density gradient, 
which is not included in that of model 1. Procedures to 
obtain the equilibrium states of hot nuclear matter by 
QMD are explained in Ref. @ for model 1 and in Ref. 
[301 ] for model 2. 

We calculate the angle-averaged radial distribution 
function §ab (j) for snapshots of QMD simulations with 
16384 (or 2048) nucleons at a proton fraction x p = 0.3 
[3l| . In calculating <7ab(j)j we place 5 3 (or 10 3 ) resulting 
simulation boxes, each containing 16384 (or 2048) nucle- 
ons, to enlarge the cutoff radius of gAB(r)', we take the 
cutoff radius to be greater than 100 fm. Then we obtain 
5y4B (?) by calculating the Fourier transform of gAB(r) 
multiplied by a Hanning window [IE , which removes rip- 
ples caused by Fourier transforming gAB(r) with a dis- 
crete change at the finite cut-off radius. 

In Fig. [3] we show (S(E V )) at a fixed temperature 
(T = 1 MeV) for QMD model 1 (left panel) and at a fixed 
density (0.0660 fm~ 3 ) for QMD model 2 (right panel), 
respectively. From these figures, we observe two impor- 
tant characteristics of (S(E V )) calculated for the results 
of the QMD simulations. One is that the peak is lower 
and broader than the results for the liquid drop models 
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TABLE I: The mass of the region where pasta phases would 
appear (M pasta ) and that where nuclei exist (M nuc ) just before 
the bounce in the typical supernova simulations for T/\ W\ = 
and 0.01. Radius means the distance of possible regions for 
pasta phases from the center of the core. Mass is measured 
by the solar mass Mq. Ratio is defined as A/ past a/A/ nuc . 



T/\W\ 
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0.8-1.6xl0 6 0.30 
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FIG. 3: (Color online) (S(E„)) at x p = 0.3 calculated for 
QMD model 1 (left panel) and model 2 (right panel). Ab- 
breviations are the same as in Fig. [Tj The left panel shows 
{S(E U )) for each pasta phase at T — 1 MeV. For comparison 
we also show (S(E„)) at 0.0165 fm~ 3 and T =1 MeV for LS- 
EOS. The right panel shows (S(E V )) for various temperatures 
at a fixed density 0.0660 fm -3 . 



shown before (see the lines for SP and LS-EOS in the left 
panel of Fig. [3]). This is due to lattice disorder and irreg- 
ularities of nuclei, which can be incorporated by QMD, 
but not by the liquid drop models. The former (latter) 
enhances the cross section of low (high) energy neutrinos, 
while they reduce that of neutrinos around peak energy. 
It is noticeable that, due to the former effect, (S(E V )) 
calculated for the QMD results has nonzero values in the 
low energy region of E v < 25 MeV in some cases while 
(S(E U )) shown in Figs. Q] and [2] vanishes in this region. 
The other important characteristic is that, as shown in 
the left panel of Fig. O the peak height does not necessar- 
ily decrease with increasing density unlike the results for 
the liquid drop models shown in Figs. Q] and [2] This non- 
monotonic behavior is observed between the phases with 
spherical nuclei and rod-like ones, and also between the 
phases with slab-like nuclei and rod-like bubbles. Let us 
first consider the latter case. In the results of the QMD 
simulations, slab-like nuclei form a layered lattice and 
rod-like bubbles form a triangular lattice. The ratio of 
the number of the smallest reciprocal lattice vectors of 
a layered lattice to that of a triangular one is 1:3. This 
larger number of the lowest order of the reciprocal lattice 
vectors for the phase with rod-like bubbles might explain 
the larger peak height compared to the phase with slab- 
like nuclei 31]. For the same reason, it is possible that 
{S(E U )} increases with a transition from the phase with 
slab-like nuclei to that with rod-like bubbles caused by 
an increase of temperature as suggested in the right panel 
of Fig. [3l Here nuclei adopt a slab- like shape at 1-2 MeV 
and nuclear structure transforms to rod-like bubbles at 
3-4 MeV, above which nuclei deform and nuclear surface 
cannot be identified. On the other hand, the reversal of 
the peak height between the phases with spherical nu- 
clei and rod-like ones would be explained by the position 
of spherical nuclei in the former being easy to fluctuate 
compared to that of rod-like nuclei in the latter, where 
the nuclear matter region is connected in one direction 
and this acts to suppress the fluctuation in the transverse 



direction. 

Finally we estimate the possible mass of the pasta 
phases in supernova cores to examine their importance 
on the core-collapse supernova explosions. We have per- 
formed two-dimensional axisymmetric simulations of ro- 
tational core collapse supernovae. The initial condition 
is taken from the profile of the 15 Mq nonrotating pro- 
genitor with solar metalicity at zero-age main sequence 
(ZAMS) by Woosley et al. [33|. The mass of the core is 
1.58 Mq. Rotation is then added in the same way as in 
Ref. [If. Here we set T/\W\ = and 0.01, where T/\W\ 
is a ratio of the rotational kinetic energy T to the grav- 
itaional binding energy W of the core. These values of 
T/|W| are typical and the latter one is also suggested 
by simulations of the evolution of massive stars with 
a typical rotational velocity on the equator at ZAMS, 
~ 200km/s [33j,|35j. The magnetic field is set to be zero. 
We use an EOS by Shen et al. [36( ■ Other basic equations 
and input physics are shown in Refs. [3~i| and [37| . 

In Tabled we show the mass M pasta of the pasta phases 
just before the bounce. The region where the matter con- 
sists of the pasta phases is defined by the following two 
conditions: 1) nuclei exist and 2) the volume fraction of 
nuclei is larger than 1/8, which is a condition for the 
fission instability of spherical nuclei giving a reasonable 
estimate of the density at which matter starts to consist 
of nonspherical nuclei [38|. Note that the EOS used in 
our simulations only incorporates the phases with spher- 
ical nuclei and uniform nuclear matter and it does not 
take account of the phase with spherical bubbles and any 
other pasta phases. Therefore values of Af pas t a in Table 
|T] are lower limits in our simulations: some region of uni- 
form nuclear matter close to a boundary with the phase 
with nuclei in our simulations corresponds to the pasta 
phases. We also show the mass M nuc of the region where 
nuclei exist in the core because the pasta phases mod- 
ify the cross section of elastic neutrino scattering, which 
is dominated by the contribution of nuclei. It is notable 
that about 10-20 % of spherical nuclei are replaced by the 
pasta nuclei in both the rotating and no-rotating cases. 
In addition, the moderate rotation in the core extend the 
region of pasta phases by centrifugal force [39(. 

In summary we have shown how pasta phases affect 
energy-dependent cross sections for coherent scattering 
of neutrinos in collapsing cores. At zero temperature, 
existence of the pasta phases instead of the phases with 
spherical neuclei and bubbles increases the energy of the 
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peak of the cross section at E v ~ 30 -40 MeV and de- 
creases the peak height. They lead to a reduction in the 
opacity of the core for low energy neutrinos of E v < 30 
MeV. At non-zero temperatures, the peak height further 
decreases but increases the cross section of neutrinos at 
energy lower and higher than that of the peak. The pasta 
phases, which amount to 10-20 % of the mass of super- 
nova cores, would modify the lepton fraction of the core. 
The next step is to construct a systematic numerical ta- 
ble or a useful fitting formula of the cross section for elas- 
tic neutrino scattering in the pasta phases and perform 
hydrodynamic core collapse simulations, which incorpo- 
rate the effects of the pasta phases on the neutrino cross 



section studied in the present article. 
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